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Abstract . 



Subroutine ’’STADET" analyzes a structure to find if 
it is stable and determinate. The structure must first be 
idealized as a rigid body supported by a system of links. 
The program checks on the arrangement and number of the 
support links. Flags are set so that this subroutine may 
be used to select whether to use simple statics or strength 
of materials methods in subsequent analysis of the problem. 
Testing of the subroutine is discussed. 

Two input-output programs are also described. One 
calculates the reactions in the equivalent support linkages 
from the resultant load; the other assumes that the pro- 
perties of the links are known and calculates the allowable 
load . 
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Notation 



A 


A matrix of coefficients 


A j B j C f 0 


The coefficients of the equation of a 
plane 


A 


A vector 


AAA 


Rectangular components of A 


B 


A column matrix 


B 


A vector 


Bx> By, B^ 


Rectangular components of B 


ll’ l2» e^c. 


Vectors made up of the direction cosines 
of the support links, link cosine vectors 


li » B j , Bj^ 


Indicates the ith, jth or kth of Bi, B , 
etc. 2 


F 


A Fortran floating point, numerical 
format 


F 


A vector force 


Fi 


The ith vector force 


Fx’ Fy, F, 


Rectangular components of F 


Fxi >Fyj_ i^zi 


Rectangular components of Fj_ 


I 


A Fortran integer, numerical format 




Unit orthogonal base vectors in the 
X, y, z directions respectively 


KPROB 


An array used to indicate the kind of 
problem involved 


1, m, n 


Direction cosines of a line 


M 


Vector moment 


Mx, My, M, 

NOSOL 


Rectangular components of M 

An array used to indicate the reason 
for the lack of a solution 
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Pl> ^2’ 


Position vectors giving the locations 
of the points of application of the 
links with respect to the first point 
of application 


PLA 


An array used to store the coefficients 
of the equations of the planes determin 
ing the lines of action of the links 


r 


The position vector of the point of 
application of a force 




Rectangular components of r 


U 


Poynting's vector 


W 


A weight 


X 


A column matrix 


X, y, Z 


Rectangular coordinates 


Xo ,7o )Zo 


Coordinates of a given point 


^1 


The scalar result of forming a dot 
product 


2i 


The vector result of forming a cross 
product 


a 


The angle included between two vectors 
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1. Introduction. 

Digital computers have been extensively applied to 
the solution of mechanics problems. However, these ap- 
plications have been limited in scope. The computer has 
been used as a "super slide rule" to solve a particular 
type of problem which was too time-consuming or too dif- 
ficult previously. It is the intention of Professor Gilles 
Cantin to attack the problem of the design of structures 
by a digital computer in a general rather than a specific 
manner . 

One of the first routines needed for the Cantin Pro- 
ject was a routine to reject unstable structures and, 
perhaps, make some selection as to the methods of analysis 
to be used in subsequent tests. The subroutine ’’STADET" 
analyzes an idealized structure to find if it is stable 
and determinate. It sets flags which can be used to print 
out a description of the kind of problem involved. V/hen 
the subroutine finds the first condition which would pre- 
clude the solution of the problem by simple statics, a 
flag is set and control returns to the input-output program. 
If the "error" is caused by a redundant support system, 
other methods of analysis may be available whereas there 
is no solution for an "unstable" problem. 

The definition of stability which is considered here 
is more restrictive than the question of stable equili- 
brium usually encountered in basic statics courses. Here 
is meant, not the criterion of energy change with small 
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perturbations, but the geometrical stability of infinitely 

strong members. As expressed by Kinney,^ 

A stable structure will support any conceivable 
system of applied loads, resisting these loads 
elastically and immediately upon their applica- 
tion the strength of all members and the capacity 
of all supports being considered infinite. 

Consider the two links supporting a weight shown in 



Fig. 1. 




Fig. 1. Support links before and after displacement. 

In Fig. la, with the links horizontal, there is no 
vertical component of the tension in the bars to support W. 
Even with infinitely strong foundations and members, the 
geometrical arrangement precludes a vertical component un- 
til some downward displacement of W has occurred (Fig. lb). 
The methods of strength of materials are required to find 
the equilibrium position of the structure. Since W is not 
resisted "immediately upon its application," this arrange- 
ment fails the test for stability. 

Both a structure and its support system must be 
stable. The problem which is examined here is limited 
to that of the external stability of a structure. The 

^J. Sterling Kinney, Indeterminate Structural Analysis, 
p. 20, 1957. 
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structure itself is considered to be a rigid body so that 
only the system of supporting links is of concern. No 
loads are imposed during the stability analysis since 
equilibrium is not involved. 

However, in some cases forces such as gravity must be 
considered in order to have a structure which meets this 
restricted definition of stability. This is a familiar 
assumption in basic statics where, for example, the gravity 
force provides a reaction against upward motion of a roller 
nest. It is up to the user of this program to decide 
whether a "gravitational” restraint is sufficient or whether 
a positive restraint is needed. 

The user must first idealize his structure so that the 
subroutine can operate on it. Appendix A, which shows the 
linkage diagrams of some common supports, may be helpful in 
this regard. However, it is at this point that engineering 
judgment is required. Needless to say, if the proper equi- 
valent support system is not employed, no reliance can be 
placed on any answers which are produced. 

Fig, 2 shows the decisions which are made by the sub" 
routine. The titles within the blocks indicate the basic 
sections of the subroutine. These titles are roughly indica- 
tive of the purpose of the sections though more information 
is usually produced than the title would lead one to believe. 

For input the subroutine requires the direction cosines 
of the support linkages, the coordinates of the points of 
application of the links and the number of links. Appendix H 
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Figo 2 Simplified Flow Chart 
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shows how the input must be arranged to be read in by the 
input-output programs of Appendix C. The output consists 
of the two flag arrays, KPROB and NOSOL. 

The subroutine first considers the points of appli- 
cation of the supporting links. If these points are along 
a line, the question of whether it is a two- or three- 
dimensional problem depends on the orientation of the links. 
If the points are in a plane, the links must all be in that 
same plane for a two-dimensional problem to be involved. If 
the points are randomly located in space, only a three- 
dimensional problem can be involved. The two-dimensional 
problem requires three linkages for stability; the three- 
dimensional problem, six. If there are fewer linkages, an 
"unstable” error print-out is made; if more, an "indeter- 
minate" one is made. 

When it has been determined that the correct number of 
support linkages is involved, the orientation of the links 
is considered. The following restrictions apply to the 
three-dimensional case: 

1 . All the links may not be parallel. 

2 . They may not lie in parallel planes. 

3. The axes of the links may not intersect one 
straight line since limited rotation would 
be possible about that line. 

4. The axes of the links may not intersect in a 
point or limited rotation is possible about 
that point. 

5. Four of the six links may not be parallel or 
intersect in any one point for, in these cases, 
a straight line can be found through the point 
(or infinity) which intersects the axes of the 

two remaining links and restriction 3 * is violated. 

In the two-dimensional case only restrictions 1 . and 



4. apply. 
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V/here the links are all parallel to a given line, 
limited motion is possible in a direction perpendicular 
to this line. 

In Fig. 3, where the pairs of links are in parallel 




Fig. 3* Arrangement with links in parallel planes. 

planes, some motion is possible in a direction perpen- 
dicular to these planes before a restoring component is 
developed. In Fig. 4 some rotation about the intersection 




■-J 

Fig. 4* Arrangement with intersecting lines of action 
of links, 

of the lines of action of the links is necessary before 
a restoring couple is developed. 

The following section and the appendices describe 
the ’’STADET” subroutine, its auxiliary subroutines and 
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the Input-Output Programs which were used to test "STADET" . 
Fortran language was used for all the programs. The testing 
was done on a Control Data Corporation 1604 digital computer. 
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2. Method. 

Subtitles indicate the sections of the flow charts 
(Figs. 5 through 9) and the "STADET'' list (Appendix J) 
under discussion. Circled numbers on the flow charts are 
used for two purposes. The first is to facilitate the 
comparison with the list of Appendix J; the second is to 
indicate a transfer to another section of the program 
when the use of a line to show the path would be cumber- 
some . 

The subroutine only operates up to the detection of 
the first error. The appropriate flag is then set and a 
return is made to the Input-Output Program. 

a. COPLANARl 

This section of the subroutine generates the 
position vectors which give the locations of the points 
of application of the links. It checks to see whether 
these points are coincident, lie on a line, are in a plane 
or are randomly located in space. In the case where the 
points of application are in a plane, this section of the 
program checks to see if the links are in the same plane 
and so determines whether a planar problem is involved. 

After deciding whether a two- or three-dimensional 
problem is involved, this section counts the number of links 
to find out if the number is sufficient for stability. 

The computer operates with a binary approximation to 
a decimal number when it is used in the floating point 
mode. Each calculation is approximate and round-off errors 



Fig. 5 Flow Chart of COPLANARl 




can accumulate., For this reason the "zero'" which must be 
used in the program is not the mathematical zero but an 
allowable tolerance. The first portion of this section 
examines the coordinates of the points of application and 
determines the value of "zero" to be used throughout 
"STADET". Appendix L discusses the method. 

Vectors are generated between the first point of 
application and the other points. These vectors are ex- 
amined to find the first non- zero one. If none of them 
is non-zero, there is only one point of support, the struc- 
ture is unstable, and a return is made to the Input-Output 
Program, 

If a non-zero vector is found, the subroutine deter- 
mines whether it is the last one or not. If it is not 
the last, the other vectors are examined to find whether 
any others are non-zero. If the first non-zero vector is 
the last of the array of position vectors, or if all after 
the first are zero, the points of application of the links 
lie on one line. 

Having determined that the points of application are 
on a line, the subroutine considers the number of links. 

A three-dimensional linkage system is unstable because the 
links intersect a line and only a planar problem has a 
unique solution. Therefore, if there are more than three 
links, the system is redundant. If there are less than 
three links, the system is unstable. If the number of 
links is not satisfactory, the appropriate flag is set. 
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V/here a second non- zero vector, P2, is found, its 
cross product with the first, P^, is formed. The cross 
product vector, is examined to find out if it is zero. 

If it is zero, P]_ and P2 are parallel. Since Pj_ and P2 
have the same initial point, this mades them coincident 
and indicates that the first three points of application 
are on a line. The subroutine then checks to see if all 
of the position vectors have been examined. If they have 
not all been examined, the search for another non-zero 
vector continues. If they all have been examined and all 
of the points of application are found to lie on one line, 
an investigation of the number of supporting links is made, 
as above. 

If a cross product vector, as found above, is not 
zero, a check is made of the number of supporting links. 

If there are only three links, there are only two position 
vectors which (since they are not parallel) determine a 
plane. If there are more than three links, there are more 
than the two vectors that were crossed to farm Xy.~ 
check whether the points of application are in a plane, 
it is necessary to dot these other vectors with X^. The 

dot product, X^, a scalar, is compared with zero. If all 
the dot products are zero, the points of application lie 
in one plane. 

With the points of application of the links in a plane 
a stable system is possible with either a two- or three- 
dimensional linkage system. If all the links lie in the 
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plane that contains the points of application of the links, 
the problem is one in two dimensions. To check for a two<- 
dimensional problem, the vectors made up of the direction 
cosines of the links are dotted with and the dot products 
are compared with zero. Since the link cosine vectors have 
one point in the point-of -application plane, a zero dot 
product indicates that the vector lies entirely in that 
plane. If a planar problem is what is involved, i.e., all 
the dot products are zero, the proper number of links is 
checked for, as before. 

When the points of application are not in a plane or 
when the links are not in the point-of-application plane, 
a three-dimensional problem is involved and six links are 
required for stability. Suitable error flags are set if 
there are not six supports. 

At this point it has been determined whether a two- 
or three-dimensional problem is involved and that the 
correct number of supporting links is involved for the 
type of problem. Next, instability due to incorrect 
orientation of the links is investigated. 

b . PARALLELl 

This section checks to see if the supporting 
links are parallel or lie in parallel planes. The first 
link cosine vector is crossed with each of the succeeding 
ones and the resulting cross product vectors are compared 
with zero. If all of the cross products are zero, all of 
the links are parallel. 
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Figo 6 Flow Chart for PAf<ALLELl 
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The second portion of this section is concerned with 
discovering whether a situation similar to the one shown 
in Fig. 3 exists or not » In the above test a non- zero 
cross product will result if two of the links are not 
parallel o (If the situation is like that of Fig. 3 > the 
cross product vector is in the direction of motion of the 
rigid body.) If all of the links are perpendicular to 
this cross product vector, the case is indeed similar to 
that of Fig. 3 » This section checks whether the links 
are perpendicular by dotting each of the link cosine 
vectors with the cross product vectors and comparing the 
resultant scalar with zero. If all of the dot products 
are zero, all of the link cosine vectors are in parallel 
planes. If all are in parallel planes, the appropriate 
flag is set; if not, the program goes on to the next 
section . 

c , PARALLEL2 

Even though all of the links are not parallel, 
the supporting linkage system is unstable, by criterion 
3 ., if four of the six links are parallel. Therefore, 
a separate check is made for this possibility. The cross 
product test is used along with a counting system. If the 
linkage system passes this test, the remaining check for a 
three-dimensional system is done by CONCURRENT!. 

d, PLANARl 

This section of the subroutine checks the three- 
link cases to see if the links are in a plane and if the 
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Figo 7 Flow Chart for PARALLEL2 
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Figo S Flow Chart of PLANARl 
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links are parallel. In Fig. 2 there are two paths from 
COPLANARl to PLANARlo The first of these three-link cases 
has the points of application on a line. In the other 
case the points of application are in a plane and it is 
knovm that the links are in the same plane. For this 
second case there is some duplication of testing, but 
the extra operating time is negligible and there was some 
convenience in prograjraning this way. 

This section crosses the first link cosine vector 
with the second c If the cross product is zero, the first 
is crossed with the third link cosine vector. If this 
cross product is also zero, the three links are parallel 
and the structure is unstable. If the second cross product 
is not zero, it is dotted with the second link cosine 
vector. If the dot product is zero, all three link cosine 
vectors are in the same plane, a stable planar problem is 
possible and the program goes to CONCURRENT!. If the dot 
product is not zero, the three links are not coplanar and 
the system is unstable. 

V'Jhen the first cross product is not zero, it is dotted 
with the third link cosine vector and the dot product is 
compared with zero. If the dot product is not zero, the 
system is unstable, as above. If the dot product is zero, 
the program goes to CONCURRENT!. 

e. CONCURRENT! 

The remaining check is to find if the lines of 
action of the links intersect at a point. In order to 
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Figo 9 Flow Chart of CONCURRENT! 
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be able to operate mathematically, it is necessary to 
convert the link cosine vectors to lines in space. The 
line form which was chosen was that of the intersection 
of two planes. (See Appendix G for discussion of sub- 
routine PLANE.) 

All of the planes are determined and the coefficients 
are stored in the array PLA. The equations in PLA are 
examined four at a time but with an overlap, i.e., equations 
one through four are examined first, then equations three 
through six; and, if a six-link problem which gives 12 
equations is involved, the process continues until all 
12 have been examined. 

The three-link system must have at least one inter- 
section or the links would all be parallel, a possibility 
which has previously been eliminated. In a six-link 
problem at least four of the lines of action of the links 
must intersect to make the linkage system unstable. There- 
fore, the system of examining pairs of lines whose equations 
are adjacent in the array PLA must find the invalidating 
intersection and a more complicated checking system is not 
required . 

The coefficients of four equations are treated as 
though they were the coefficients of four simultaneous 
equations in four unknowns. They undergo the matrix 
singularity test of REACl (See Appendix F). If the de- 
terminant of the matrix is not zero, the set of equations 
is consistent and the first two lines do not intersect. 
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If the next overlapping set of four is also consistent, 
three lines do not intersect and the. link support system 
is satisfactory., This would complete the checks by "STADET". 

On the other hand, if the determinant is zero, the set 
of equations is dependent or inconsistent and the two lines 
intersect or are parallel. If the lines intersect, the 
point of intersection can be found by solving three of the 
equations simultaneously. If the lines are parallel, the 
four combinations of these four equations, taken three at 
a time, are all inconsistent. No account is taken of this 
parallelism since it is already known that less than four 
links are parallel. The next overlapping set of four equa- 
tions is then considered. 

When an intersection is found, the coordinates are 
substituted into the equations of all the planes. The 
number of equations which are satisfied is counted. In 
the case of a planar problem, if six are satisfied, all 
the lines of action intersect and the linkage system is 
unstable. In the three-dimensional case, if eight or 
more are satisfied, the linkage system is unstable. V/here 
instability exists, the appropriate flag is set. 

This completes the checks made by "STADET”. The 
rigid body is supported by a system of links whose re- 
actions can be determined by simple statics or the appro- 
priate flag has been set in NOSOL. 

The one remaining possible error, that of not having 
the load in the correct plane for a planar problem, is not 
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strictly of concern in deciding whether the system of 
links is stable and determinate « The check for this 
error is discussed in Appendix C, Input-Output Program 



3c Testing of Program. 

In order to devise a logical system of tests the sim- 
plified flow chart of Fig. 2 was used. This shows the 
questions that the various sections of the subroutine 
answer. In all, there are 2$ distinct paths that can be 
followed. Only four of these paths lead to a satisfactory 
problem, i.e., one which can be solved by simple statics. 
Fig. 10 lists the types of problems which were solved. The 
complete print-outs of the solved problems are collected 
in Appendix K. 

Problems 26 through 29 are the complements of the 
satisfactory problems of the first 25. The ”bar reactions 
output of TESTIA was used as the "allowable bar loads" of 
TESTIB. The output of TESTIB, the "allowable load torsor-t 
should be, and is, the "resultant torsor" of TESTIA. 

In problems 11 and 17> where there are more than six 
links, the storage process is such that the last link is 
stored on top of the first. This results in a nonsense 
print-out of the structure data. In problem 11, for ex- 
ample, the first link has two direction cosines of one. 

In the use of TESTIA to solve a planar problem it 
is necessary to consider the possibility of having an 
invalid problem because the load is not in the same plane 
as the structure. Since this possibility is caused by an 

■'Torsor -- A coined word (combination of torque and 
tensor) used to indicate a generalized resultant which 
is composed of a couple and a force. 



Fig. 10 
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error on the part of the one setting up the problem and 
IS not strictly a factor in the stability of the structure, 
testing for this possibility was made part of the Input- 
Output Program o The testing of a planar problem therefore 
requires one step beyond "STADET". The subcases of prob- 
lems 22 and 23 show the results of such tests. 

It was not possible to get an accurate estimate of 
the time required to solve a problem. It takes roughly 
one second to solve an unsatisfactory problem and about 
twice as long to solve a satisfactory one. This differ- 
ence is not in the computation time but in the time re- 
quired to write the output on the magnetic tape. In fact, 
when a series of problems is being solved, only a slight 
pause can be detected between the end of writing one 
answer and the beginning of the output for the next. To 
these computing and writing times must be added the time 
to compile the program. This is about 75 seconds for 
TESTIA and about 15 seconds for TESTIB. 
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4o Conclusion o 

When this project was undertaken, it was understood 
that this subroutine would eventually be used as part of 
a much more ambitious program., The objective of providing 
a routine which decides whether the structure may be 
analyzed by the methods of simple statics has been achieved. 

For the purposes of this routine a decision was made 
to have storage for only six equivalent links, the maximum 
number compatible with a solution obtainable from statics. 
This is the only restriction which must be removed in 
future, expanded use. The arrays KPROB (kind of problem) 
and NOSOL (no solution) were set up to be used as keys to 
select the proper methods of attack » If NOSOLfl) is set 
to one, a redundant structure is involved and more power-= 
ful methods must be used to effect a solution. KPROB 
indicates whether the simplification of a planar problem 
is valid o 

While it is probable that a bona fide planar problem 
would be set up in one of the coordinate planes, it might 
be advisable to provide a subroutine to transform the 
coordinate system so that an oblique plane problem is 
transferred to a coordinate plane. It may or may not be 
desirable to transform the results back to the original 
coordinate system. 

The next step in the generalized analysis of structures 
by digital computers could be along one of several paths. 
Since this subroutine is concerned wlch only the external 
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stability of the structure, a subroutine to examine the 
internal stability needs to be written » The simplifica- 
tion of considering only space frames might yield a good 
first approximation u 

Another path would be to either use a code or have 
an examining type of subroutine to decide the type of 
problem and then jump to one of a set of subroutines 
designed to solve a given class of problems c 

It is felt that this is one small step in removing 
some of the drudgery from engineering design. 
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APPENDIX A 



Linkage Diagrams of Common Supports 



A link is a member capable of supporting only an 
axial force. It has spherical hinges at both ends so 
that it cannot transmit any moments. In the drawings 
below the link is shown thus: , where the dot indi- 

cates the hinge at the point of application. The other 
end of the link is to be understood to have a hinge which, 
in turn, is attached to a frame of reference. 

As discussed in Section 1, Introduction, some of 
these supports are not the equivalent of a complete link. 
In these cases it is the user's responsibility to either 
modify the support or to make sure that the reaction is 
in the allowed direction. 



Smooth plane support 




Because friction is neglected, 
a linkage system which only 
opposes vertical motion is 
obtained. The link shown may 
only have a compressive reaction. 



Flexible cord 




The reaction is along the cord. 
Only a tensile reaction is 
allowable . 



2S 



Ball and socket 




A ball and socket is the 
equivalent of three links 
acting at a point . 



Built-in support 




BODY 



A built-in support is the 
equivalent of six links. 

The transparent body is 
shown supported by one of 
many stable and determinate 
arrangements of six links. 
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APPENDIX B 



Mathematical Tests 



The examination of the system of supporting links is 
accomplished by the application of a few simple mathemati- 
cal tests to the vectors which are part of or are generated 
from the input data, 

Ac The dot or inner product of two vectors A and B, 
neither of which is zero, is zero if and only if A and B 
are perpendicular. 

A • B = |A| |B| cos a where a = included angle 



A * ® ^x®x ^y^y ■^z^z 



subscripts indicate x, 
y, z components of A 
and B 



Bo The cross or outer product of two vectors A and 
B, neither of which is zero, is zero if and only if A and 
B are parallel or collinear. 

A ^ B = |A| |B| sin a u where u is a unit vector per-= 

pendicular to both A and 
B and with them forming 
a right-handed set 



A X B = 



i i A 

^x ^y ^z 
By 



i, 2 A unit, orth- 
ogonal base vectors in the 
X, y, z directions, re- 
spectively 



Co The question of whether of not lines intersect is 
solved by examining the equations of the planes determining 
the lines. If a set of four equations determining a pair 
of lines is inconsistent or dependent, the lines are parallel 
or intersect. The ’’matrix singularity test” of REACl (see 
Appendix F) is used to find if the equations are dependent 
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or inconsistent » 

If the lines intersect, the point of intersection 
can be found by solving three of the four equations simul- 
taneously. 
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APPENDIX C 

Input“Output Programs 



These control programs provide for reading in data, 
operating the required subroutines and printing results. 

In both cases the first step is reading in the data on 
the equivalent linkages. Then ’’STADET" is called and the 
stability and determinateness of the support system is 
examined. The rest of the data for that problem is then 
read in regardless of the results of "STADET”. The 
intermediate operations of the two programs will be con- 
sidered separately. 

TESTIA 

If the structure is stable and determinate, the re» 
sultant of the forces and couples acting on the structure 
is calculated. This establishes the B vector in the matrix 
equation, Ax = B. The A matrix, that array which, when 
used to preraultiply the link reaction vector, x, gives B, 
is found next. If there are six links, A is six by six; 
the six equations represented are independent and the values 
of the link reactions are found by subroutine REACl. 

In the case of a planar problem, "STADET’’ has examined 
the structure and found it satisfactory, but it remains to 
be seen whether the load is applied in the correct plane. 
Therefore, checks are made to see that the resultant force 
vector is in the plane of the structure and that the re- 
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sultant couple vector is perpendicular to that plane. 

If there are only three links, A is six by three; 
only three independent equations exist among the six, 
and it is necessary to solve various sets of three equa- 
tions until a consistent set is found (See Appendix F, 
REACl). When a consistent set is found, the results are 
printed out as discussed below. 



TESTIB 

If the structure is stable and determinate, the A 
matrix is calculated. The allowable loads in the various 
links are knowns in this problem. Therefore, the x vector 
is multiplied by A giving the allowable load on the struc- 
ture. In this case there is no question about the plane 
of the resultant load as in TESTIA. 

Output 

In both programs the intermediate calculations are 
suppressed if the structure is not stable and determinate. 
The printout then consists of the input structure, a 
description of the problem type and the error. If the 
structure is satisfactory, the input for the Intermediate 
calculations and their results are printed. Since the 
case of coincident support points would be an unusual 
blunder, a special error printout was not used. Since 
"STADET” is not really operated in this case, the print- 
outs are: 
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NOT A PLANAR PROBLEM and 

STRUCTURE UNSTABLE LINES OF ACTION OF SUPPORTS 
INTERSECT A LINE 

(See problem 1.) 

Both programs reset conditions so as to solve the 
next problem and then continue until all solutions have 
been made. (See Appendix H for problem format.) 
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PROGRAM TESTIA 

DIMENSION VB(6,3),B(6,3) , FORCE! 100,3) , CORO (1 GO, 3 ) , COUPL ( 50 , 3 ) , 
1KPRCB(5) , NOSOL (6) ,RFCR(3) , RCCPL ( 3 ) , BR EAC ( 6 ) , CCMX (6 , 7 ) ,PS(6,7), 
2ARI3) ,BR(3),CRCSS t 3) , VP ( 5 , 3 ) , PO INT ( 3 , 4 ) , PLA ( 1 2 ,4 ) , PL ( 4 ) , PM ( 4 ) 
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51COFORMAT(5X, 61HSTRUCTURE UNSTABLE LINES OF ACTION OF SUPPORTS ARE 
lONCURRENT) 
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35 F0RMAT(3(3X,F1U.7) ) 
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309 IF (Jl-5) 310,311,311 

310 J1=J1+1 
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IF (NCO) 314,315,314 
314 CONTINUE 



PRINT 125 

PRINT 135, ( (COUPL ( I , J J J ) , JJ J = 1 , 3 ) , 1=1 ,NCO) 
PRINT 91 
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B( I, J ) 
571 VB(I,J) 
DO 572 



PRINT 136 

IF(NP-NPROB) 520,580,580 
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PROGRAM TESTIB 

DIMENSION V8(6,3) ,B(6,3) , KPROBI 5 ) , NOSOL ( 6) ♦ RFOR ( 3) ,RCOPL(3) 
1BREAC(6) ,COMX(6,7 ) , AR ( 3 ) ♦ BR ( 3 ) , CROSS ( 3 ) , VP( 5 » 3 ) , PO INT(3,4 ) , 
2PLA( 12,4) ,PL(4) ,P M ( 4 ) , BFORCE I 6 ) ,PS(6,7) 

COMMON VB,B,NBAR, KPROB, NOSOL tCOMX , VP, BFORCE 
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NOSOL( I )+l 
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IF (KPROBM)) 157,158,157 
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APPENDIX D 



Subroutine CPOST 



The subroutine evaluates the cross product of two 
vectors expressed in rectangular coordinates. A Fortran 
expression is written for each term of the expansion of 
CROSS. 



CROSS = 



1 i Is 



\ 

Bx By Bg 



= (AyB^-A^Byli + (A^Bj^-A^Bg) j + (A^By-AyB^)k 
= CROSS(l)i + CR0SS(2)j + CR0SS(3)k 
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SUBROUTINE CROST ( AR , BR , CROS S ) 
DIMENSION AR(3) ,3R{3),CROSS(3) 
CROSSm = AR(2)»BR(3)-AR{3) »BR(2 ) 
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APPENDIX E 
Subroutine RESl 



The subroutine finds the resultant of a system of 
forces and/or couples by first adding up the rectangular 
components of either and then calculating the effect, if 
any, of the forces on the resultant couple. 



F = 


£ 




etc 
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XI 
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F . 

XI 


etc 



r - + Tyj + r^k 



F = F^l + Fyi + F^k 

M = r X F 

= + Myi + M^k 



is the X component of F^ 
is the X component of 



r is the radius vector to the 
point of application of the 
force 



r^, r„ and r„ are the rectan- 

JL y ^ 

gular components of r 

F is the vector force 

F^, Fy and F^ are the rectan- 
gular components of F 

M is the vector moment 



M^, My and Mg are the rectan- 
gular components of M 
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SUBROUTINE RESI 

DIMENSION VB(6, 3) ,B( 6»3) , FORCE! 100, 3) , CORD ( 100,3), COUPLI 50,3) , 
1KPR03 (5) , NOSOL (6) , RFOR ( 3 ) , RCOPL { 3 ) , BREAC { 6 ) , COMX ( 6 , 7 ) , PS ( 6 ,7 ) , 
2ARI3) ,BR(3) ,CROSS (3) ,VP(5,3) , POI NT { 3, 4 ) , PL A ( 1 2 , 4 ) , PL ( 4 ) , PM ( 4 ) 
COMMON VB,B,FORCE , CORD, COUP L, NFO,NCO,NB AR,KPROB, NO SOL, RFOR, RCOPL 
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70 CONTINUE 

IF(NFO) 80,90,80 
80 CONTINUE 

DO 30 1=1, NFO 

RCOPK 1 )=RCOPL { 1 ) -FORCE! 1,2) *CORD( 1,3) +FORCEI 1 , 3)»CCRD( 1,2) 
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APPENDIX F 



Subroutine REACl 

This subroutine is a modification of the GAUSS2 sub- 

2 

routine used by C. B. Bailey in his program for solving 
simultaneous linear equations. Whereas Bailey was solving 
the matrix equation, Ax = B, with a maximum of 50 vectors 
B for each coefficient matrix A, nowhere in this program 
is more than one B vector involved with a given A. 

The routine locates the largest first column coef- 
ficient and, if necessary, exchanges the row containing 
this coefficient with the first so the largest element is 
. Multiples of row one are subtracted from the other 
rows to make the column one coefficients below the first 
zero. Next the coliimn two coefficients below the first 
are examined to find the largest and rows are exchanged, 
if necessary, as before. This Gaussian elimination process 
is continued until all the elements below the main diagonal 
are zero. 

At each reduction step, the value of the diagonal 
element is compared with zero. A very small diagonal 
element indicates ill-conditioned equations and causes 
an error output. This is the ’’matrix singularity" test 
used in CONCURRENT!. If the matrix is non-singular, a 
back- solving process is used to find the unknowns. 

^C. B. Bailey, F2 UTEX LINEQN, CO OP Manual, Dec., I960. 
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DIMENSION AA(6,7) ,X(6) 
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IF (AA(K,D) 32,3U,32 



32 RATIO = AA(K,L) /AA(L,L) 
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APPENDIX G 



Subroutine PLANE 

Given a vector with direction cosines 1, m and n 
and one point, the problem is to find a pair of planes 
whose intersection is parallel to the vector and passes 
through the point . 

First, the subroutine checks to see if the vector 
is ’'nearly" parallel to a coordinate axis. "Nearly’’ is 
defined as having a direction cosine greater than O.G, 
i»e., lying within about 37° of the coordinate directions. 
If the vector is "nearly’’ parallel to one axis, it is 
crossed with vectors along the perpendicular axes. 

If the vector is not "nearly" parallel, a check is 
made to see whether it lies within about 6° of a coordinat 
plane, i,e., the direction cosine is less than 0.1, If 
it is close to a coordinate plane, one of the crossing 
vectors is chosen perpendicular to that plane. 

If neither of the above tests is met, the given 
vector is crossed with vectors in the x and z directions. 

There are now two cross product vectors of the form: 

= A_l + B j + Ck 

Since these vectors are perpendicular to the planes con- 
taining the vectors which formed the cross products, the 
dot product of X^_ arbitrary vector in one of the 

planes is zero. An arbitrary vector is of the form: 
(x-Xo)i + (y-yo)i + (z-z„)k 
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I 



where (x^,y^,z^) is the known point. The dot product has 
the form: 

A(x--x^) + B(y-=y^) + C(z-z^) = 0 
This may be rearranged to: 

Ax + By + Cz = D 
where: D = Ax^ + By^ + CZq 

The coefficients A, B, C, and D are stored in the array 
PLA. 



6l 



I 



i ■■ 

i. -I 






'■j 

.■i 






,v 4 flWcrt/ ^1'^* * '(,•*• -.X » d*x®iiw'' 

, fflrtoljiad.? 

r. i.^ , , - ' ,<-icr 

! / -. ‘.-a T «>cr 'C'l^ aidlT 



. ’ r 

,/,r *^X3 4a‘19dV 



.LVi 



SUBROUTINE PLANE! CC,DD,PL,PM) 

DIMENSION CC( 3) ,00(3 ) ,PL ( 4) , PM( 4) , AA( 3) , BB( 3) 



ro 



O' 



m 

M 

Q 



CO 

I 



ro ‘— 















ro 










































f~ 


%— 




























r— 


o 






o 




• 


• 




• 


• 




o 


u 


o 


o 


o 


U 




o 




II 


o 










r— 


•— 
















II 


X— 


UJ 




»— 


UJ 














II 


II 


il 


II 


II 


II 




LL 


:o 




a. 


Z) 


II 


II 


J- 


II 


II 


-3- 














:»£C 








lo 


2T 




























CO 




. — 


00 








o 






o 










-cr 




*— 


< 


»— 


r~ 


< 


k- 


•— 


(O 


k- 


r— 


CNJ 


k— 






W' 


w> 


w> 






w> 






w> 


?• 


X—^ 








wo 




< 


-J 


•C 


CO 






C 


LL 


c:) 


o 


u. 


O 


< 


cc 


o 


< 


CO 


o 


< 


Q. 


Cl 


CO 


CL 




Q 




u 


o 




O 


< 


CO 


o 


< 


CO 


o 






i 












r~- 






r- 


uO 






o 







C\i 



62 



AA( 2) 



5 



w 

O w 

o o 

jc< a 







. — . 


















-j 




















a. 


CL 




-J 
















#> 




Cl 


2! 












< 


OD 






CL 












< 


cc 




+ 


+ 












•» 


•» 










#• 


•» 






o 


o 


ro 


. — . 














o 


o 


•» 


-n- 


St 




o 


un 










f— 








•» 


•> 








h- 




«J 


jr 




to 


CM 






oO 


on 


M 


a. 


CL 




•ft 


•ft 




LU 


O 


O 










<M 


o 






a: 


q: 


::cr 


II 


II 




w 








o 


o 








2: 














o 


• — . 




or 


O 


c 






-j 


-j 


»— 




sr 


3 


h- 


I— 






«j 


-j 






wr 


»— 






Q 


a 


< 


< 


o 


-j 




LU 


o 


o 




u 


o 


o 


o 


CL 


CL 


or 


o 


o 


UJ 


j* 






00 




O 




rn 


O' 





63 



APPENDIX H 



Problem Format 



TESTIA 



Fortran Format 


Symbol(s) 


15 


NPROB 


Hollerith 




15 


NEAR 


3F10.a,3F14.7 


VB,B 


13,12 


NF0,NC0 


6F12.6 


FORCE, CORD 


3F14.7 


COUPE 



Input data 

Total number of problems 
Problem title = up to ^0 
characters 
Number of links 
Direction cosines, coordi- 
nates of links 
Number of forces, couples 
Force components, points of 
application 
Couple components 



Repeats from the Hollerith title for each problem « 
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TESTIB 



Fortran Format 


Symbol ( s ) 


Input Data 


15 


NPROB 


Total number of problems 


Hollerith 




Problem title - up to SO 






characters 


15 


NBAR 


Number of links 


3F10.S,3F14.7 


VB,B 


Direction cosines, coordi= 






nates of links 


6F12.6 


BFORCE 


Allowable link loads 


Repeats from the 


Hollerith 


title for each problem 
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APPENDIX I 



STADET Flags 



KPROB = Kind of Problem 

1 Support points are on a line 

2 Support points are in a plane 

3 Reaction matrix is singular 

4 Planar problem (if 4 is 1) 

Not a planar problem (if 4 is 0) 

5 Resultant load is not in plane of structure 



NOSOL 

1 

2 

3 

4. 

5 

6 



No solution because: 

Structure unstable - lines of action of supports 
are concurrent 

Structure unstable - lines of action of supports 
are parallel 

Structure unstable - insufficient number of 
supports 

Structure indeterminate - excess number of supports 
Structure unstable - lines of action of supports 
intersect a line 

Structure unstable - supports are in parallel 
planes 
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APPENDIX J 



STADET List 
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'j..' TvlCAtr 




SUBROUTINE STADET 

DIMENSION VB(6,3) , B ( 6 , 3 ) , FORCE (1 00 , 3 ) , CORD (1 00, 3 ) , COUPU 50 , 3 ) , 
1KPR08(5) , NOSOL (6) , RFOR ( 3 ) , RCOPL ( 3 ) , BRE AC ( 6 ) , COMX ( 6 , 7 ) , PS ( 6, 7 ) , 
2AR(3) ,BR(3) .CROSS (3) ,VP( 5,3), POINT (3,4) ,PLA (1 2, 4 ) , PL ( 4 ) , PM ( 4 ) 
COMMON VB, B, FORCE .CORD, COUPL, NFO.NCO, NBAR , KPROB , NO SOL , RFOR, RCOPL 
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IF (I-NBAR ) 23,315,23 
21 IF (I-NBAR+1) 24,25,25 



26 IF {ABSF(VP(K, 1 ))+AeSF(VP(K,2) )+ABSF( VPIK.3) )-EP) 32,32,31 

32 IFIK4-1-NRAR) 33,25,25 

33 K = K<-1 
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53 IF ( J-K)54,52,54 

54 C =CRCSS( 1 ) *VP{J, n + CR0SS(2 ) *VP( J , 2)+CR0SS( 3 )«VP(J ,3) 
J1 = Jl+1 
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END 



APPENDIX K 



Problem Types 

In this appendix "bars” and ''supports" are used as 
synonyms for "links". See Section 3» Testing of Program, 
for a discussion of the problems. 
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PLANAR PROBLEM - TOO FEW SUPPORTS 
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.REDUNDANT STRUCTURE - MORE THAN THREE BARS 
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UNSTABLE THREE BAR STRUCTURE - BARS ARE PARALLEL 
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UNSTABLE THREE BAR STRUCTURE - BARS ARE NOT IN A PLANE 
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STRUCTURE UNSTABLE INSUFFICIENT NUMBER OF SUPPORTS 
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STRUCTURE UNSTABLE LINES OF ACTION OF SUPPORTS ARE CONCURRENT 



REDUNDANT STRUCTURE - MORE THAN THREE BARS 



oooooo 
o o o o o o 
oooooo 
oooooo 
oooooo 
oooooo 
oooooo 

M • • • • • • 

oooo 

) i 



o 



< 

0 oooooo 
^ oooooo 
-J oooooo 
a oooooo 

01 oooooo 
< oooooo 

oooooo 



o 



00 



o 

a. 



UJ 


oooooo 

oooooo 

oooooo 

oooooo 




cc 


oooooo 




o 


oooooo 






oooooo 




o 


X • « • « • • 




o 


ooooo 




a: 


r-f— CNJ— eg 




H- 

00 

O 

I 

s 

o 


oooooo 






oooooo 




-j 


oooooo 




o 


00 oooooo 




LL 


0:2:000000 






< oooooo 


UJ 


UJ 


CG oooooo 


3 


X 




< 


h- 


LL 


-J 




0 




OF 


t >0 


< 


> 


uJ 

3 oooooo 






oooooo 


»— 'X 


►—* 


uo oooooo 


UJ 


-J 


0 oooooo 


LLI-J 


-H 


0X000000 


CiiU 


CO 


oooooo 


<0 


< 


3 oooooo 


DC 


h- 


CJ ••••♦• 


UOCL 


00 


»—4 


h- 




h- 


3c< 


UJ 


0 


>«.< 


X 


UJ 


03 


H- 


q: 






oooooo 


«J 


UJ 


Cl oooooo 


V-CL 


3 


oooooo 


c: 




oooooo 


o< 




^oooooo 


. CL 


c 


oooooo 


CLH- 


X 


oooooo 


co 


UJ 







STRUCTURE INDETERMINATE EXCESS NUMBER OF SUPPORTS 



UNSTABLE STRUCTURE - SUPPORTS ARE PARALLEL 
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UNSTABLE STRUCTURE - LINES OF ACTION OF SUPPORTS INTERSECT 
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STRUCTURE UNSTABLE LINES OF ACTION OF SUPPORTS ARE CONCURRENT 



10. UNSTABLE STRUCTURE 
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11. REDUNDANT STRUCTURE - MORE THAN SIX BARS 
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12. UNSTABLE STRUCTURE - SUPPORTS ARE PARALLEL 
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STRUCTURE UNSTABLE LINES OF ACTION OF SUPPORTS ARE PARALLEL 



13. UNSTABLE STRUCTURE - SUPPORTS ARE IN PARALLEL PLANES 
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14. UNSTABLE STRUCTURE - FOUR OF SIX SUPPORTS ARE PARALLEL 
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15. UNSTABLE STRUCTURE - LINES OF ACTION OF SUPPORTS INTERSECT 
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17. RGDUNDANT STRUCTURE - MORE THAN SIX BARS 
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18. UNSTABLE STRUCTURE - SUPPORTS ARE PARALLEL 
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19. UNSTABLE STRUCTUIIE - SUPPORTS ARE IN PARALLEL PLANES 
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20. UNSTABLE STRUCTURE - FOUR OF SIX SUPPORTS ARE PARALLEL 
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21. UNSTABLE STRUCTURE - LINES OF ACTION OF SUPPORTS INTERSECT 
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22a. PLANAR PROBLEM 
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22b. PLANAR STRUCTURE l^ITH IMPROPER LCAO 
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23a. OBLIQUE PLANE PROBLEM RESULTANT FORCE IS IN PLANE 
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23b. OBLIQUE PLANE PROBLEM RESULTANT FORCE IS NOT IN PLANE 
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24. STABLE SPATIAL PROBLEM 
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THE BAR REACTIONS ARE 

39.0014903 82.3859488 44.5109861 -28.6600889 -39.1484822 -47.350468 



25 . STABLE SPAT I AL . PROBLEM 



oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

ooorooo 






o 

o 

o 

o 

o 

M • 



>or^ 



2: 

o 



< 

o 



oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

tnooooo 

soOLDOi^r- 



a. 
cl 
< 

LL>- 

O ^ CT'G" 



00 



I- O 
<r o 
o 

Q rv) 
LU *— 
»— • >- • 
-J C\i 

a 

a 

< 






00 

00 

00 

00 

oco 

oro 











0 






CNJO 




0 






0 






• • 




a 






0 






00 










0 






r-o 




oooooo 






h- 






a- 




oooooo 






X • 






1 




oooooo 






a 








UJ 


oooooo 




00 


1 








cc 


oooooo 




0 








W'W* 


3 


omcMOOo 




< 




0 






h- 


oot-ootn 




0 




0 




4- 4- 


0 


X « • • • • • 




a 




0 






3 










0 




33 


cc 


— CM 




0 




0 






K- 






X 




rvjO 




■ — — 


00 








0 






00 








3: 


0 


• 




00 


0 






0 


0 


0 




00 




*. 




a 


0 


f— 




00 








a 


CsJ 






00 


3 T 






0 


M • 






0 — 


0 


r-.— — 000 




a 


ao 






r-o 


-J 


^^^000 






•— 






• • 


3 


oooooc> 




a 










0 


00 CMCMCMOOO 




X 








m— 


a 


0 ^ Xjj- J* j-ooo 




a 








0 




< (NJCNJCNJOOO 








0 




— 


UJ 


CO .d* J* O'*— 




a 




0 






X 






a 




0 






a 


a 1 1 








00 0 








0 




(V' 


0 


a 0 






a 






0 


0 


a>-Lr) 


*» 


4. 4. 


a 


00 




00 


0 


as:— 


00 






UJ 




cc 


00 0 


3 


1 -^ 


►— M-H 


> 


2 000—00 




0 


a r-* 


a 0 






h- 


— 0004 T 00 


2 i. 


a 


o> • 


0 0 


cc 







00 oooooo 


a 




aao 


O' 


0 


00 


3 


0 OOOOJOO 


a 


a 


0 fO 




00 


00 




02 i. 000 .d -00 


a 


ZL 


a 




cc 


00 


GO 


OOOCSJOO 


0 


< 






0 


00 


< 


^ — OOj-'-O' 


cc 


a 






a 


00 




0 •••••• 


a 


a 








CNjrvo 


00 


1 — « 




3 






a 


— CO 




a 


Qi 


00 




0 


z 


• • 


a 


0 


< 


a 




0 


< 


0— 


X 


a 


'Z. 


a 


0 


0 


a 


fO 


»— 


CC 


< 




0 


0 


a 






— 0000— — 


a 


a 


0 


0 


3 




iXJ 


0 COOO.d^ 


a 


X 


<\j 


xo 


CO 




X 


oooooo 




a 


f— 


xo 


a 




1 — « 


OOOOCNJCsl 


< 




X • 


• 


cc 




5 : 


aooood'.d 




Q 


ao 








< 


OOOOCvKM 


a 


Z 


f >0 




a 


H li 


X 


0^— — 0 ^^.d 


0 


I— • 






X 




UJ 




X 


a 






a 


aX 



106 



THE BAR REACTIONS ARE 

-121.0548607 -106.1962450 -32.8508693 41.9365735 89.0738227 167.042582 



26 . PLANAR PROBLEM 
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27 . OBLIQUE PLANE PROBLEM RESULTANT FORCE IS IN PLANE 
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29. STABLE SPATIAL PRCELE.'^ 
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APPENDIX L 



Definition of ’’Zero” 

The computer must convert any non-integer number to 
an approximate binary form which will fit its word length. 
The conversion and length limitation allow an accuracy of 
about ten decimal digits. With each calculation round-off 
errors accumulate. Whereas the flow charts ask if a 
quantity is zero, what must really be considered is 
whether the quantity is less than the maximum expected 
round-off error. If there are no conflicting requirements, 
the maximum round-off error may be estimated and, with 
perhaps some leeway, defined as "zero". However, where 
different orders of magnitude may be involved in the same 
problem or in subsequent problems or where it is advisable 
to have a large tolerance for one calculation and a small 
tolerance for another, the situation becomes far more 
complicated. One cannot define the "zeros" in terms of 
the number of significant figures in an expression; a 
specific magnitude is required. 

An attempt was made to tie the size of the "zero" to 
the size of the structure. It was felt that barring un- 
reasonable complications, such as having a very small 
structure at a great distance from the origin^ the system 
of averaging the coordinates of the points of application 
and multiplying this average by the computer accuracy 
would give a suitable zero. However, instead of being 
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able to use 10 times the average coordinate, it was 
necessary to use a "zero" 10,000 times as big to get 
proper answers. Further investigation revealed that 
the difficulty probably involved the use of subroutine 
REACl in CONCURRENT!. 

In subroutine REACl a small value of "zero" allows 
small diagonal elements which, during the back-solving 
process, generate large coordinates of the point of 
intersection. When these large values are then used 
in CONCURRENT!, the round-off error becomes much greater 
than the "zero" defined above. Unsuccessful efforts were 
made to modify the "zero" obtained above in a logical way 

for use in CONCURRENT! and REACl. 

-6 

A multiplier of 10* gives proper answers, at least 
for this set of problems, but there is no other justifi- 
cation for its use. The author, regretfully, did not 
have sufficient time to pursue this investigation. 
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